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Abstract Various aspects of the magnetism of the quiet sun are reviewed. The 
suggestion that a small scale dynamo acting at granular scales generates what 
we call the quiet sun fields is studied in some detail. Although dynamo action 
has been proved numerically, it is argued that current simulations are still far 
from achieving the complexity that might be present on the Sun. We based this 
statement not so much on the low magnetic Reynolds numbers used in the simu- 
lations but, above all, in the smallness of the kinetic Reynolds numbers employed 
by them. It is argued that the low magnetic Prandtl number at the solar surface 
may pose unexpected problems for the identification of the observed internetwork 
fields with dynamo action at granular scales. Some form of turbulent dynamo at 
bigger (and deeper) scales is favored. The comparison between the internetwork 
fields observed by Hinode and the magnetism inferred from Hanle measurements 
are converging towards a similar description. They are both described as randomly 
oriented, largely transverse fields in the several hecto-Gauss range. These similar- 
ities are ever making more natural to assume that they are the same. However, 
and because of the large voids of magnetic flux observed in the spatial distribution 
of the internetwork fields, it is argued that they are not likely to be generated 
by dynamo action in the intergranular lanes. It is concluded that if a dynamo 
is acting at granular scales, the end product might have not been observed yet 
at current spatial resolutions and sensitivities with the Zeeman effect. Thus an 
effort to increase these resolutions and polarimetric sensitivities must be made. 
New ground- and space-based telescopes are needed. The opportunity offered by 
the Solar Orbiter mission to observe the Quiet Sun dynamics at the poles is seen 
as one of the most important tests for confirming the existence, or otherwise, of a 
granular ly driven surface dynamo. 
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1 Introduction 

A consensus about the existence of a small-scale dynamo (SSeQ) operating at 
the solar photosphere is being consolidated in today's solar physics (see, e.g. , 
IVoeler and Schusslerll2007l : lAbbettl2007t IPietarila Graham et al ]2010HSteiiJ2012L 
for a review). High resolution magnetograms from ground and (mostly) space- 
based telescopes observed in t he internetwork are often used to in dicate that such 
a surface dynamo exist (e.g. IPanilovic et al 1 l2010al lLitedl201l|) . On theoretical 
grounds, simulations of various kinds have been use d to suggest how universal 
the various ingredients of such a dynamo seems to be (|Moll et al.ll201ll) . They all 
indicate that turbulent shear stresses acting on the inertial range act as the main 
mechanism able to efficiently convert kinetic into magnetic energy. The conclusions 
from the various simulations of turbulence in a conducting fluid seems to be that 
it would have been a lot harder to explain the absence of a surface SSD than its 
presence. With an emphasis on the observational side, we review in this work the 
current status of this consensus and try to pinpoint which aspects are more solidly 
established and which are less settled. 

Section [2] will be the only one that concentrates on simulations and the theo- 
retical aspects related to the problem of the existence of SSDs on the Sun. It will 
address the most controversial argument questioning the existence of such a mech- 
anism, namely the fact that the solar convective zone has a magnetic Prandtl num- 
ber that is orders of magnitude smaller than one, while simulations work near the 
~ 1 regime most of the time. Low magnetic Prandtl numbers are known to severely 
discourage dynamo action since the early simulations of convectively driven tur- 



bulent dynamos (jNordlund et al.lll992tlSchekochihin et al.l l2004a. 2005). Recently, 



progress has been achieved, however, that indicates that an SSD is indeed possible 
in the low magneti c Prandtl regime (jlskakov et al.ll2007t ISchekochihin et alJl2007t 
lBrandenburdl201ll) . But the situation is not conclusive and the papers addressing 
this issue often resort to the fact that a mixed polarity field is observed at the 
solar surface as the firmest indication that such a mechanism should exists. How- 
ever, and in the absence of a clear proof that this observed (internetwork) fields 
originate from an SSD -and such a proof is not available yet— the only progress to 
settle this issue will come from further work in the simulation front. 

The observational arguments that have been put forward to favor the pres- 
ence of a solar surface dynamo are discussed later in Section [3] There are ba- 
sically two such arguments. First, the evidence from the observed Hanle effect 
in lines such as the Sri 4607 A line and the careful modeling of these signals 
indi cate the existence of a ta ngled field with a mean strength of < B >~ 130 



G ([Truiillo Bueno et al.M2004j) some few hundred kilometers above the solar sur- 



face (see also Truiillo Bueno et al.l 120061 and references therein) . This number 



was originally derived under some model assum ptions that made it uncertain 
to w ithin a factor two. However, a recent study (jShchukina and Truiillo Buenol 



2011 ) of the predicted Hanle signals from the MHD simulations described in 



Pietarila Graham et al.l ((2009b) has eliminated some of this model dependency 



and confirmed such large mean field strengths. As a mixed polarity tangled field 
at unresolved scales leaves basically no trace in the Zeeman profiles, these fields 



1 Small scale here refers to generation of magnetic fields at scales smaller than the energy 
injection one, the granulation. 
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have always been a prime candidate to be considered as originated from a surface 
SSD. If the tangling occurs at scales near or above present resolutions, some signa- 
tures can be detected, tho ugh. It is unclear if these fields have been observed using 
the Zeeman effect (but see lLites et al.ll2009l : lBellot Rubio and Orozco Suaredl2012[ 
and Section[3|). The second observational argum ent in favor of an SSD comes from 
the H inode spectropolarimeter (SP) instrument ( Tsuneta et al.l2 008: Ko sugi et al.l 
120071) and its unprecedented characterization of the internetwork fields using the 
Zeem an effect (|Lites et al.l20 08; Ishik awa and Tsunetal201lllOrozco Suarez and Bellot Rubiol 
l2012h . While transverse fields were known to exi st in the quiet su n, as originally 
found by the Advanced Stokes Polarimeter fASP: lLites et al"1ll996l) . it was totally 
unexpected that these fields have a predominant transverse character. This trans- 
verse nature seems to fit in a natural way with an origin related to a turbulent 
dynamo as shown by recent simulations (see ISchussler and Voglerl[2008h . 

A possible outcome given this state of affairs could be as follows. The existence 
of an SSD acting at granular scales at the solar surface can eventually be confirmed 
from a set of improved SSDs simulations (along the lines described in Section [2]). 
In them, a continuous distribution of fields is obtained that is able to explain the 
Hanle depolarization levels from those fields created at the smaller scales and the 
largely horizontal internetwork fields from those at larger scales. The separation 
between these two sets of fields does not have to be sharp and a range of spatial 
scales can contribute to both the Zeeman and Hanle results (or, perhaps, that 
the two observed processes are due to fields exactly at the same scales as it can 
be inferred from the recent results of lBellot Rubio and Orozco S uarez 2012]). Note 
that, such a field distribution would solely depend on the existence of the always 
present turbulent convective motions near the surface and, thus, should be inde- 
pendent of latitude and of activity cycle phase. While this conclusion seems rather 
plausible given the current evidence, the aim of this work is to address some of 
the known problems that might prevent such an outcome. In particular, Section 
B] describes some observations recently obtained with the IMaX/SUNRISE mag- 
netograph (|Martmez Pillet et alJkoiltlSolanki et al.ll2010f) that show regions that 
display very little magnetic activity, either measured as residual signals in time- 
averaged deep magnetograms or as evidenced by a lack of flux emergence episode s 
in the form of small-scale loops (as discovered by iMartmez Gonzalez et al1l2012l) . 
It is unclear how these voids are compatible with a granular ly driven SSD. That 
the situation is far from clear has been corroborated recently by the study of 
(|Stenfloll20ll based on SDO/HMI magnetograms) who proposes the existence of 
a basal flux of order 3 G that is suggested to be an upper limit to the efficiency 
of an SSD at the solar surface. According to this result all the internetwork fields 
observed with Hinode/SP will not be generated through such a mechanism and 
only the Hanle depolarizing fields could be originated through it (if at all). 

In spite of this somewhat confusing situation, it is important to stress that 
our understanding about the nature and the properties of the quiet sun fields has 
improved enormously in recent years. But it is clear that a number of important 
questions remains on the theoretical/modelling side and on the observational front. 
Section [5] finishes this work proposing a way forward to further advance in this 
understanding of the quiet sun magnetism. Not surprisingly, we promote an effort 
to increase the polarimetric sensitivity and the spatial and temporal resolutions of 
both, the Hanle and the Zeeman observations. Studying the statistics of the quiet 
sun fields at various latitudes will also prove crucial. 
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2 Small scale dynamo action at low P m . Implications for the solar case 

The seminal reference that triggered the present debate on the existenc e of a con- 



JTs 



vectively driven turbulent dynamo at the solar surface was the work of Cattaneo 
1999), although the debate is older (see, e.g., IPetrovav and Sza kaly 1993; Lin 



1995). It is important to point out that this work mentioned, both, the granular 
and supergranular scales as possible contributors to such non-helical dynamo. The 
simulations used closed upper and lower boundaries with vertical fields in both 
of them. The Reynolds and magnetic Reynolds numbers could be clearly defined 
for this simulation thanks to the fixed computational grid used to solve the MHD 
equations. They were Re = ul/v = 200 and Re m = ul/r] = 1000, with / the char- 
acteristic length of the energy injecting convective cells, u the velocity of these 
cells and v and r/ the molecular viscosity and the magnetic diffusivity, respec- 
tively. These numbers are large enough to ensure the development of turbulence. 
But it is important to note that Re was five times smaller than Re m (magnetic 
Prandtl number of P m = Re m /Re = vjr\ = 5 > 1). Under these circumstances, 
the magnetic field sees a smooth mean flow efficiently acting on it. The numerical 
simulation resulted in dynamo action saturating at 20% of the kinetic energy flow. 
The crucial ingredient was the chaotic nature of the driving flows. Figure 2 of this 
paper already showed that, at the surface, the strong fields were localized in the 
downflow lanes, while the cell interiors showed no (vertical) field signature. The 
situation was different in deeper layers where fluctuating fields were filling basi- 
cally the whole volume. At the time of the publication, the dominant transverse 
nature of the internetwork was not known (see Section [3]) and this aspect was 



not an alyzed. For this reason, the profiles synthesized by Sanchez Almeida et al.l 



(2003) using these simulations concentrated on the study of the asymmetries in- 



duced in the Stokes V profiles (circular polarization) and its comparison with those 
observed in the internetwork. A shortage of asymmetries indicated that the sim- 
ulations still did not achieve as much complexity as present in the Sun. However, 
using this synthesis, and after including effects such as telescope diffraction, it was 
predicted that when improving the spatial resolution from 1 arcsec to 0.15 arcsec, 
one should detect four times more Stokes V polarization signals. 

A number of assumptions made in this simulation (such as the Boussinesq ap- 
proximation) have been relaxed in more recent works. The more realistic ones (in 
terms of their proximity to the physical conditions o n the Sun) are those made 
with the MURaM code (IVorier and Schiisslerl l2007t ISchiissler and Voglerl 120081; 
iPietarila Graham et al.l |2010T) . In particular, thev have addressed the important 
questi o n of the role played by the closed boundary conditions assumed in lCattaneol 



1999). IStein et al.l (|2003h pointed out that, in the solar convective zone, fields are 
submerged efficiently to the bottom of the convective zone by strong and concen- 
trated downflows showing little recirculation near the surfac e. This recirculation 



was artificially enhan ced in the simulatio ns of lCattaneol ( 1999() by the use of closed 



boundary conditions. IStein et al.1 ((2003) concluded that diverging upflows sweep 
the fluid into downflows, often vortical, where stretching and twisting becomes ef- 
fective (and balanced by diffusion) but all these fields are rather rapidly submerged 
down into the bulk of the convective zone. The energy added to the flux that vis- 
its the surface was a very small fraction of the global budget of magnetic energy 
and the effect cannot be considered a local dynamo. This criticism has, however, 
been superseded by the MUraM simulations which used an open boundary and 
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allow for a non-zero pointing flux at the bottom boundary. The way in which this 
boundary condition was implemented in the simulations of Voaler and Schiissler 
(|2007t ) was by imposing an artificially increased magnetic diffusivity there. This 
diffusivity ensured that horizontal fields moving downward in the simulation leave 
the box unimpeded while, at the same time, prevented horizontal flux from en- 
tering the domain. All flux leaving the bo ttom boundary in these simu lations was 
created by dynamo action inside the box. IVogler and Schusslerl (|2007|) concluded 
that while the downward pumping of flux outside of the domain does indeed reduce 
the growth rate of the dynamo, it does not shut it down. They predict that as long 
as a sufficiently high Re m is used (above the critical magnetic Reynolds number, 
Re^), dynamo action will be observed in the simulations. Specifically, they see 
exponential growth of the magnetic energy for Re m of 2600, while they find a 
decrease for Re m below 1300. Now, here we should caution that the exact values 
of Re m and Re (and thus of P m ) are not as easily defined as in ICattane o (1999) 
simulations. The need to resolve shocks with artificial viscosity schemes or the 
implementation of the bottom boundary diffusivity necessarily implies that these 
numbers are more difficult to ascertain. An estimate of th e effective Reynolds and 
Prandtl numbers for the MURaM simulations was given bv lPietarila Graham et all 
(|2010[ ) using various moments of the velocity and magnetic spectra (or Taylor mi- 
croscales). The various dynamo runs available from this code turn out to have 
Re TO 6 [2100,8300] (with the latter value using a grid resolution of 4 km) and 
Pm £ [0.8,2]. This set of simulations all displayed an SSD generating magnetic 
fields inside their volume. The conclusion was that current MURaM simulations 
(with P,„ ~ 1) will show dynamo action as long as Re m > Re m ~ 2000. The main 
mechanism identified for this generation was the stretching and twisting of field 
lines by fluid motions in the inertial range of the spectrum of velocity fluctuations. 
In particular, it was concluded that dynamo action is concentrated in the turbulent 
dow nflows as fie ld line stretching against magnetic tension is very efficient there 
(e.e. IStem1l2012|. for a review). In order to clarify the nature of the observed dy- 



Moll et all ( 201 lh have analyzed the underlying physical mechanism 



namo action, 

under various physical conditions and assumptions. They concluded that for the 
cases studied (incompressible MHD, Boussinesq convection and compressible solar 
convection), the field is amplified by similar inertial range shear stresses that are 
independent of the conditions at the injection scale. The inclusion of compress- 
ibility effects or the asymmetry between upflows and downflows generated by the 
strong stratification did not influence the result. They, thus, termed the dynamo 
mechanism as universal. 

One concern remains about the existence of a possible SSD on the Sun, though. 
The problem has been known since the early studies of dyna mo action in con- 
ducting fluids. It was originally formulated by iBatchelorl ( 1950l) who studied how 
turbulent motions stretch the field lines and amplify the magnetic energy as long 
as this process remains unimpeded by ohmic diffusion. Field line stretching is pro- 
duced by fluid motions in the inertial range whose dissipation scale is set by the 
viscosity of the fluid. High viscosity v couples the field lines to the flow and allows 
it to bend them efficiently. Large magnetic diffusivity 77 decouples the plasma (and 
the flows) from the field lines and prevents the bending. Thus, it was always clear 
that the efficiency of the SSD was going to be controlled by the interplay of these 
two effects as measured by the ratio P TO = v/i]- Clearly, both conditions Re m 2> 1 
(to favor field line stretching) and P m 3> 1 (to couple fluid motions and field lines) 
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Fig. 1 The Reynolds number, magnetic Reynolds number and (i nverse) Prand tl number in 
the solar convective zone according to the mixing-length model of Spruit (1973)- The X near 
the surface marks the tipical values achieved of both Re and Re m in the simulations. 



boost local dyn a mo ac tion. Based on an analogy between vorticity and magnetic 
fields iBatchelorl f|l950h even concluded that for P m < 1 no SSD was possible. This 
conclusion was later criticized by a number of authors as the analogy cannot in- 
clud e the different initial and bo undary conditions seen by these two fields (see, 
e.g. Boldvrev and Cattaneoll200 4h . However, simulations in the early 90's (see, e.g., 
iNordlund et alJll992j) . including compressibility and strong stratification, already 
resulted in dynamo action only for P m > 1 with efficient shutting down of the 
dynamo for P m < 1. Other simu l ations encoun tering the same problem are dis- 
cussed in lBoldvrev and Cattane o ( 2004) and in ISchekochihin et al.l ( 2005f ). Thus, 
the question of the existence of an SSD at low Prandtl numbers has received some 
attention in recent years. 

Before we briefly describe the results from the numerical studies about the 
existence of an SSD at low Prandtl numbers, it is important to remember what 
are the actual numbers that occur in the solar convective zone and get an idea of 
how far or how close are we from simulating these conditions. To this end, Figure[T] 
shows Re, Re m and P m as computed in the mixing-length based model o f lSpruitT 
(jl974l) . The velocities at the injection scale u are obtained as part of the model 
and the characteristic length is assumed here to be I = z/2, with z the depth 
inside the convective zone. The magnetic Reynolds number changes from 10 5 at 
the photosphere to 10 9 at the bottom of the convective zone, whereas the Reynolds 
number stays constant at a level of 10 12 . This makes P m g [10 -7 , 10~ 3 ], with the 
smallest value reached at the photosphere. Thus, everywhere in the Sun, we have 
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Re 2> Re m 3> 1 and P m <C 1 which is exactly the regime where the existence of 
an SSD becomes problematic. It is generally believed that for a sufficiently high 
Re m there will always be an SSD at work. But the above mentioned simulations 
prompted a deeper study about the existence and nature of an SSD under solar 
conditions. As already mentioned, the simulations bv lCatta neo (1999j) had P m = 5 



and those from the MURaM code always move close to the P rn ~ 1 case (there is 
one with P m = 0.8 that is discussed below). Note that the regime where Re m and 
Re are similar is actually very favorable for the numerical codes as similar grid 
sizes resolve the dissipative scales of both, magnetic and velocity fields. 

What is the physical argument behind this difficulty to generate an SSD when 
Pm <C 1?- Under such conditions the viscous {l v ) and resistive (Iri) scales follow 
ln/l v ~ P m 3/4 > 1 (|Schekochihin et alJl2~004al. l2005t) and the resistive scale l n falls 
in the middle of the inertial range. Turbulent eddies of scales I > do the nec- 
essary field line bending and twisting for dynamo action at a rate of ui/l (with 
u\ the typical flow velocity at this scale). If P m > 1 only these eddies occur and 
the field lines always see a spatially smooth flow. However, if P m < 1 one has 
now eddies below the resistive scale I < l v . These eddies act on the field as a 
turbulent diffusion with diffusivity u;Z and destroy magnetic energy. It is the pre- 
dominance of this last process what can make dynamo action impossible as it was 
seen in the previously mentioned si mulations. Using an incom pressible spectral 
MHD code and the PENCIL code 0, ISchekochihin et all (|2005l) studied what are 
the possible asymptotic limits when Re S> Re m and the corresponding values of 
the Re m for the existence of a dynamo. The above descri bed effect always trans- 
lates i nto a sharp increase in Re m as P m — > (see also IPietarila Graham et al.l 
l2009al) . but does not prevent the existence of a dynamo in this regime. The two 
asymptotic limits are, first, as Re — > oo, Re m — > const, so that dynamo action is 
possible for higher Re m and, second, Re m — ► oo with Re m /Re — > P m =const, in 
which case no dynamo is possible (turbulent diffusion efficiently dissipates mag- 
netic energy at small scales) . Which exactly of the tw o asymptotic limits prevails 
has been under much debate in recent years. While ISchekochihin et aD ( 2004al) 



and lSchekochihin et al.l (|2005j) (see their Figure 2) favor the existence of a P m (no 
dynamo) from simul ations of incompressible magne toconvection reaching values of 
P m as small as 0.15. lBoldvrev and Cattaneol (2004) provided analytical arguments 



favoring the existence of a Re m . 

While the debate in the mid last decade did not look promising for confirming 
the existence of an SSD at the solar surface, the s i tuatio n has ch anged recently 
(even if not completely settled). ISchekochihin et al.l (|2007{) (see also llskakov et all 
2007) performed simulations of incompressible MHD turbulence reaching values of 
Pm = 0.1 and with a sufficiently high Re m that indicated a plateau region where 
a Re m is observed (see their Figure lb). Admittedly, the number of such simula- 
tions proving the existence of this plateau is very small but the autho rs consider 
it enough numerical certainty. The results from iBrandenburgl ( 201ll) using the 



PENCIL code (that includes compressibl e effects) and low values of P m resulted 
in dynamo action being activated as well. IPietarila Graham et ail (|2010|) also find 



growth of magnetic energy in the one case they analyzed with P m = 0.8. Thus, 
the most advanced existing numerical simulations of small-scale dynamo action in 



See |http:// www. nordita. dk / software /pencil-code 
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turbulent MHD currently favor the occurrence of such a process in the low P m 
regime. However, a number of caveats remain: 

— First, and most importantly, Re m increases with de creasing P m . The exact 
factor depends on the specificities of the simulations. iBoldvrev and Cattaneol 
( 2004j) suggest a factor 7 increase in Re m when shifting from the P m 3> 1 



to the P m <C 1 case. They use an analytical model of isotropic and homoge- 
neous tur bulence that includ e s the extra roughness of the velocity field for 
P m < 1. ISchekochihin et al.l (2007) find from incompressible forced turbu- 



lence a factor 3 increase. As in the MURaM simulations one has Re m ~ 2000 
( Pietarila Graham et al.ll2010|) . this means that as soon as we move into the 



low-P m regime, we need magnetic Reynolds numbers above at least 6000 
to be able to trigger dynamo action. These magnetic Reynolds numbers are 
not currently a c hieved by this code. In particular, none of the runs used by 
iDanilovic et al. I (|2010al) would be able to actually sustain dynamo action. On 



top of that, the saturated field strength is known to decrease with decreasing 
P m , thus the expected field strengths will be smaller than those computed for 
P m ~ 1. The exact amount of this reduction is still a v ery controversial issue 
(see ISchekochihin et al.ll2004al 120071: lBrandenbur"gll201ll) and its magnitude for 
the solar case unknown. But the net effect will be a reduction in the fields as 
compared to those computed for P m > 1. 

To complete the demonstration of the existence of a dynamo driven by fluid mo- 
tions in the inertial range at low P m values, a growth rate of the dynamo scaling 
with R e^ 2 must be obtained from the sim ulations. Neither [Schekochihin et all 
(2007j) nor IPietarila Graham et all ( 2010h have reached that (see Figure 3 of 



the latter work). Simulations with an increa sed resolution are needed to finally 
settle this issue. ISchekochihin et "ah | (|2007l) concludes that, as long as this is 



not achieved, the mechanism that sustains the growth of the magnetic field 
fluctuations in the l ow-P m re gime will rema i n bas ically unknown. 



fluctuations m tire low-f m regime will remain basically unknown. 
S cheko c hihin et al. | (|2007|) and llskakov et al.l (|2007|) concluded that in the P m <S 



1 regime, the magnetic energy spectra is fundamentally different from that 
found in the P m S> 1. The spatial distribu tion of the growing magnetic fields is 
qualitatively different too (see Figure 2 in ISchekochihin et al.ll2007h . This indi- 
cates that the use of simulations in the P m ~ 1 range to compute the ensuing 
Stokes profiles and its comparison with those observed may not be justified. 
The often simulated case with P m ~ 1 has th e same spectral energy pro perties 
and field distribution as the P m S> 1 case ( Schekochihin et al.l 2004b). This 



i s pro bably why the case with P rn = 0.8 simulated by IPietarila Graham et all 



i s pro bably wny tne case witn f m = U.o simulated by|rietar 
(2010) was so similar to the those in the range of P m 6 [1,2] 

Let us finalize this section by stressing that if the situation looks confusing, it is 
because this has indeed been th e case in this topic for some time (see llskakov et all 
l2007t ISchekochihin et aTll2007l. who speak about a frustrating outcome). One ar- 
gument commonly given to promote the existence of an SSD is the observation 
of a mixed polarity field in the internetwork regions of the Sun (and many of the 
above me n tioned wor ks use this argument one way or another). Figure 2 of, both, 
ICattaneol (|l999|) and IVoeler and Schusslerl (|2007h clearly indicate that this is a 
reasonable argume nt. The point we want to s tress here is that the same applies 
to Figure 10 from IStein and Nordlundl (|2006l) . which does not include an SSD. 



IPietarila Graham et al. ( 2009a|) estimate for this simulation a Re m ~ 600, which 
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is known to be too small to develop dynamo action. However, their mixed polarity 
distribution lo cated in the interganular li nes looks as 'solar' as in the other cases. 
In the work of IStein and Nordlun d ( 2006) emphasis is made on diverging upfiows 
bringing flux to the surface, expulsion to the intergranular lanes and sweeping 
of field lines into strong downflows that carry the flux into deeper layers. These 
simulations extent typically further down than those that concentrate in SSD 
generation and also include larger scales such as those associated with the meso- 
granulation. As shown in Fig. [1] the deeper we move into the Sun, the larger Re m 
and P m (although still smaller than one). Thus, a valid question is if it is not more 
natural to ask if a solar SSD exists at meso- and supergranular scales and, if so, 
whether they dominate over that that might exists at granular ones. This point 
will be further discussed in Section [3] 



3 Observed signatures of small scale dynamo action at the solar surface 

We now turn to the observational aspect of the discussion and ask the question: 
Have we seen the fields produced by a possible SSD operating at the solar surface? 
As it will become evident, there has been, as before, solid observational progress 
and areas with much confusion. Basically, two candidates exist that are often 
considered as by-products of an SSD, the internetwork fields and the, so-called, 
hidden fields that generate the Hanle depolarization signatures. 



3.1 Zeeman signals 

There is no question that Hinode/SP data has produced a major quantitative 
and qualitative jump forward in o ur understanding of the internetwork fields. The 
publication by iLites et al.1 (|2008l) of a gigantic slit scanned map with consistent 
10 polarimetric sensitivity and homogeneous spatial resolution of 0.3 arcsec 
changed our view of the quiet sun magnetism. In this map, a myriad of patches 
with predominant linear polarization signatures was discovered. The ASP already 
found the existence of episo dic burst of larg ely transverse fields (the Horizontal 
Internetwork Features, HIF. lLites et al.lll996l) but they were thought to be rather 
sporadic. The only previous indication of their existence and g lobal character came 
from the SOLIS instrument as found bv lHarvev et al.l (|2007t) . But no prediction 
from the simulations or estimate of their magnitude was available. The existence of 
this ubiquitous ho rizontal field has now r eceived full confirmation from the SUN- 



PJSE/IMaX data (|Danilovic et al.ll2010bf) who cou ld make the first study of their 



evolution (see the animation in Solanki et aljfioioh and establish a solid statistics 



of their lifetimes. Both instruments, Hinode/SP and SUNRISE/IMaX coincide in 
locating these HI F at the borders of the upflowing granules for a large fraction of 
their evolution. ILites et~ I (|2008|) emphasized that the linear polarization signa- 



tures were not co-spatial with line-of-sight fields that were more frequently found 

in the intergranular lanes. 

Before the publication of lLites et al.1 ( 2008f) most of the discussio n on the inter- 
network flux concentrated on obtaining its mean unsigned flux fsee ISolanki| [2009. 
for a review on this topic previous to the impact of the Hinode measurements) . 
The idea was that an intricately complex field with mixed polarities observed 
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in the best available Stokes V magnetograms was the outstanding description of 
the internetwork fields. Increased resolution (or sensitivity) will result into ever 
increasing amounts of longitudinal signals observed as there was less and less 
cancellation due to instrumental effects. Either because of lack of reliable mea- 
surements of the field inclination or because of a habit to focus studies of solar 
magnetism exclusively in longitudinal magnetograms, no mention to its possible 
transverse character was traditionally made. As the field stren gths were expecte d 
to be near or below equipartition values (less than 500 G, see lKeller et al.lll994j) . 
these fields were not thought to be necessarily vertical either. Before the Hinode 
results, measurements of < \B^\ > (the spatially averaged unsigned longitudinal 
flux in the internetwork) were rou tinely being made and its v ariation with the 
spatial resolution closely followed (|Sanchez Almeida et al.l l2003h . In a way, Hin- 
ode/SP results have made this emphasis obsolete. We now know that these fields 
are largely transverse and one should mainly care about < \Bt\ > or simply about 
the spatially averaged < |Bj >. We should caution here that these magnitudes are 
obtained by observations of different Stokes parameters that have different sensi- 
tivities to the real magnetic field components on the Sun and to the fraction of the 
observed pixel that they occupy (the filling factor). Thus, the steps to compute 

< \B\ > from the observed < \Bj\ > and < \Bt\ > are more problematic than 
what one might anticipate. iLites et al.1 ( 20081 ) estimated that the quiet sun map 
obtained by Hinode/SP had a < \B app \ > of 11 G (or Mx cm" 2 ) and a < \B T PP \ > 
of 55 G. These estimates were based on using integrals of the Stokes parameters 
that were calibrated against magnetic fluxes but with no account for the fraction of 
the pixel occupied by the fields. This is why they are named 'apparent' fluxes. An 
analysis performed bv lOrozco Suarez et aD ((2007) of the same data, but this time 
using a Milne-Eddington (M-E) inversion code including a filling factor as a free 
parameter, resulted also in a predominantly transverse nature of the internetwork, 
albeit with a smaller ratio of transverse to longitudinal appare nt fluxes. 

It is no exaggeration to say that the ratio measured by ILites et al.l (2008) 

< \B T PP \ > I < \B'^ PP \ >~ 5 came as a surprise and was, thus, subjected to a deep 
scrutiny by the community. 

Several factors can create a systematic bias in this ratio. Spatially averaging 
noise affects a positively defined quantity such as Bx in a way different than 
what it does to a signed quantity (Bl). As the Zeeman effect has a sensitivity 
different for each of these two components, the visibility of a given field strength 
is different depending on whether it is a field aligned with the LOS or perpen- 
dicular to it. In particular, fields close to the noise levels translate into different 
visibility thresholds. Last but not least, there is the already mentioned difference 
in how the filling factor couples to the real field strengths and inclinations for 
each one of the two components and depending on the specific method of analysis 
used. These effects and their impact into the factor 5 obta ined by the first Hin- 
ode/SP measurements have been studied by various authors (lAsensio Ramodl2009l: 
Borrero and Kobell l201ll IStenflol boilt ISanchez Almeida and Martinez Gonzalez! 
20 lit ISteiner and Rezaeill2012L see the latter for a review). It all translates into 
understanding how exactly noise influences the final result given the method of 
analysis one follows and the various thresholds for inclusion of a given pixel or 
not. Depending on the specific case, different values for < \B T PP \ >, < >> 
or inverted parameters B, inclination, azimuth and filling factor ar e obtained. 
Figure [2] shows the central portion of the same magnetogram used in ILites et al.l 
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Fig. 2 Continuum intensity from the central portion of the iLites et al.l II2008T ) Hinode/SP 
map (top left), signed magnetogram from the same portion (top right), unsigned magnetogram 
(bottom left) and transverse magnetogram (bottom right). 



(2008) with the top panels displaying the continuum intensity and scaled 
to ±20G. Network bright points are visible (e.g., top-right of the Figure) and the 
corresponding large Stokes V signals evident in the magnetogram. The bottom 
two images provided IB^ 13 ] (left) and \B^F P \ (right) both scaled within the same 
range [0,140]G. Almost all of the signals seen in the left image are identified with 
network regions as can be recognized by inspection with the two top panels. Thus, 
they will be excluded when computing < \B^ PP \ > for internetwork regions. In 
contrast, all of the signals seen in the right panel correspond to internetwork and 
contribute to the < |-By PP | > average. But note also in this last panel that regions 
with no apparent signals do not show up as dark black as in the unsigned longitu- 
dinal image, but they show a grey shade. This background is created by noise and 
it is larger than for the longitudinal magnetogram simply because of the effects 
described above. These regions should neither be just included when computing 
< \B^F P \ > nor completely excluded (as some transverse fields might exist there). 
It is in all of these details that reside some of the contradictory numbers that 
have been published. The most promising venue to clarify the situation is that 
of reducing noise. With the current instrumentation this can only be done by us- 
ing torigjsrjjitejrr^ has b een r ecently pursued very successfully 
bv | Otqzco Suarez and Bellot Rubia ( 20121) and iBellot Rubio and Orozco Suared 
(|2012h . They use slit integrations of 6.1 minutes that reach a polarimetric sensi- 
tivity close to 10 -4 , one order of magnitude better than commonly achieved. For 
visible lines, this improvement lowers down the detectability threshold for average 
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transverse fields by a factor 2-3. lBellot Rubio and Orozco show that 

slit positions with almost 10 minute effective integration time that harbor linear 
polarization signals basically everywhere (60% above 4.5<r). 

The increased exposure times impose a penalty in the sense that the spatial 
resolution is decreased and evolutionary effects are intermingled in the final results. 

This should not give the impression that these signals are always present on the 
Sun and detectable if sufficient sensitivity is available as they might have occurred 
in these pixels for only a fraction of the exposing time. Nevertheless, the data ob- 
tained with these long exposures is perfectly suited to reduced the noise induced 
bia s present in the previous analysis of Hi node/SP . The M-E inversions performed 
by Orozco Suarez and Bellot Rubiol (|2012|) and iBellot Rubio and Orozco Suarezl 
(|2012h . consistently show that the internetwork fields have intrinsic field strength 
typically in the range of 100-200 G with basically no kG present, a strongly peaked 
field inclination distribution near 90°, with most of the pixels displaying inclina- 
tions in the range [45°, 135°], and an azimuth with no preferred orientation. Filling 
factors move in the range of [0.2,0.4]. The authors inverted only those pixels that 
had a sufficiently large linear polarization signal (in either Stokes Q and/or U) to 
ensure a reliable result from the inversion. These results are largely free from most 
of the concerns expressed so far on the nature of internetwork fields. One criticism 
that remains to this analysis is the use of the M-E approximation and inversions 
able to reproduce the asymmetries are desirable. But the main conclusions from 
these recent analysis are likely to be confirmed by these more complex inversions, 
as those based in the M-E approximation are known to pr ovide robust atmospheric 
mean s even in the presence of complex stratifications fsee lWestendorp Plaza et al.l 
1998). Thus internetwork fields do not have an isotr opic distribution of inclination 
as it has been argued in a number of recent works ( Asensio Ramosl [20091 IStenfld 
l201ll : ISanchez Almeida and Martinez Go nzalez 2011 ) and their field strengths are 
typically on the few hG range. 

Inversion codes allow an inference of the filling factor as a se parate free parame- 

ter (a dmittedly, the most model dependent of all of them). Thus lOrozco Suarez and Bellot Rubiol 
(2012) were able to give real mean (not apparent) fluxes. For the real fluxes, they 
obtained, < \B T \ >~ 198G, < \B L \ >~ 64G and < \B\ >~ 220G (here the <> aver- 
age means those pixels with large enough Stokes Q and U signals to allow a proper 
inversion, not the complete map). Their ratio, now, is < \B^\ > / < \Bl\ >~ 3.1. 
Note that because this ratio has eliminated filling factor effects and the inclusion 
of only those poi nts that were inve rted, it is not directly comparable to the num- 
ber provided by iLites et al.l (|2008h . However, it confirms the largely transverse 
nature of the internetwork fields as originally shown in that work. By eliminating 
filling factor effects, this new ratio allows for a cleaner comparison with numerical 
simulations. 

Interestingly, although not anticipated, numerical simulations of magnetocon- 
vection seem to have no problem in generating large amounts of horizontal fields. 
Soon after the publication o f the results from Hinode, the S SDs simulations from 
ISchiissler and Voglerl (|2008h (MURaM code) and those from lSteiner et al.l (|2008j) 
explained that these large amounts of transverse fields were present in thei r sim- 
ulation boxes at different heights. The simulations from lSteiner et all ( 20081) were 
not dynamo simulations and used instead imposed fields in both vertical and hori- 
zontal directions. These two initial conditions generated a predominant horizontal 
field in the region where the Fel line pair observed by Hinode forms. Thus, while 
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SSDs are capable of generating a domin antly transverse fie ld, it is not an exclu- 
sive property of them. In the work of ISteiner et al.1 (|2008|) . it was through the 
well known flux expulsion mechanism of vertical fields to the intergranular lanes 
that horizontal field lines were expelled above the granules in the overshooting 
region, generating the predominant transverse fields. The quantitative comparison 
with the observed fie lds was more complicat ed and both simulations fell short of 
the values observed. iDanilovic et al.l ( 2010al) used the MURaM simulations and 
performed spectral synthesis including instrument degradation and n oise to com- 
pare the values predicted from the simulations with those observed by ( Lites et al.l 
|2008|) . The result was that while the factor five in the ratio of apparent mean fluxes 
came naturally out of the SSD simulations, the absolute flux levels were close to 
those observed only if the SSD fields were multiplied artificially be a factor 2-3. 
After this artificial increase, the average values in the SSD simulation are closer 
to 100 G in the formation region of the Fel lines which is n i cely c ompatible with 
the peak field strength in Orozco Su arez and Bellot Rubiol (|2012|) . Similarly, the 
simulations presented bv lSteiner et alh 20081) resulted in average fields of the order 
of 20 G for the transverse component and suffer from the same problem as the 
SSD simulation. One could argue that the small Re m numbers achieved in the SSD 
simulations and the field strength introduced in the simulations were too low and 
simply increasing them will explain the higher fluxes encountered by the observa- 
tions. In any case, what was clearly established from all these studies was the fact 
that the natural state of a magnetic field component below equipartition strengths 
and closely coupled with the solar granulation is that of a predominant transverse 
field component as found by Hinode/SP. Another conclusion is that while SSDs 
are compatible with this result, the latter cannot be offered as a demonstration of 
their existence at the solar surface as a non-dynamo magnetoconvective simulation 
found the same results. 

Another attempt to investigate whether the internetwor k fields are generated 
by dynamo action has been presented recently bv lLited (|201ll) . In this work, the po- 
larity imbalance of the internetwork field regions in 45 Hinode/SP maps is studied. 
Arguably, polari ty balance i s considered a necessary outcome of an SSD. However, 
and as stated bv lLited ( 201ll) . this is not a sufficient condition as after a sufficiently 
large number of turnover times, the same polarity balance is to be expected in non- 
dynamo simulations. The difficulty here stems from the fact that to measure the 
internetwork polarity imbalance one needs to carefully isolate these fields from 
network ones. If the internetwork fields originate somehow by the shredding of 
nearby network fields, one expects the internetwork to have the same sign of the 
imbalance and a proportionality between the two. Interestingly enough, while no 
scaling with the unsigned flux was found, a suggestive correlation between the 
signed flux imbalance in the internetwork and that in the nearby network was 
measured (see Figure 3 of the paper). However, being conclusive with this result is 
difficult as the isolation between internetwork and network fields is always prob- 
lematic. This approach deserves further study, probably including some stray-light 
correction that deconvolves the wings of the spread function. This would allow to 
decontaminate the internetwork fields from the surrounding network contribution 
and allow a more reliable study of the resulting polarity imbalances. 

It is clear that it will be very difficult to conclusively demonstrate that the hG, 
predominantly transverse internetwork fields originate from a granulation driven 
SSD. As mentioned in the above paragraph, we can only aim at disproving the 
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SSD hypothesis rather than expect a firm confirmation of its presence. In Section 
4, we present a recent result that, if consolidated, could be considered as one such 
refutation. 



3.2 Hanle signals 

Scattering polarization in spectral lines and its modification via the Hanle effect 
allows to study a completely different parameter spa ce of solar magnetism not 
accessible with the Zeeman effect (see the reviews in iTruiillo Bueno et al.|[200l 
Ide Wiin et~aLll2009t fStenfktolll) . A tangled field distribution within the resolu- 
tion element is invisible through the Zeeman effect but can leave a clear imprint 
in the linear polarization profiles as long as the field strengths are below the Hanle 
saturation values (typically, a few hundred Gauss). While strong homogeneous 
vertical fields are hardly hidden in Stokes V through the Zeeman effect, weak 
disorganized transverse fields show up easily in Stokes Q and U of selected spec- 
tral lines thanks t o the Hanle effect. T his different sensitivity of Hanle effect has 
provided recently (|Truiillo Buenoll201 1) a seemingly alternative description of the 
Quiet Sun fields to that di scussed above (and inferred from the Zeeman effect). 
Ever since the early studies (|StenfldfT 982'l. the Hanle signals have been interpreted 



as being due to a tangled mixed polarity field. For this reason, it has always been 
natural to associate these fields with the outcome of a tu rbulent dynamo (see, e.g. 
IVogler and Sc hiissler 2 0071: IPietarila Graham et al.lfeoiol) This alternative descrip- 
tion can be summarized as follows. The center-to-limb variation of the depolariza- 
tion signals observed in the optically thick Sri 4607 A line suggests the presence of 
a tangled field with characteristics strengths of B ~< 130 > G (as mentioned in the 
Introduction). This value is computed using realistic atmos pheric models of the 
solar surface and complex 3D radiative tr ansfer calculations! Truiillo Bueno et al.l 
120041 : IShchukina and Truiillo Buenoll201l[) . In contrast, the analysis of the depo- 
larization signals of a set of optically thin molecular lines suggests a mean field of 
the order of only < 10 > G (|Truiillo Bueno et alJlioollKleint et al.ll201ll ). These 
two distinct fields can be made compatible if one realizes that the molecular lines 
are entirely formed in the hot smooth upflowing granules an d that these lines are 
blind to the fields present in the turbulent intergranular lanes (jTruiillo Bueno et al.l 
l2004fh The Sri line however sees the fields in both regions and, in order to give 
rise to average values of 130 G over the whole solar surface, ITruiillo Bueno et al.l 



( 20041) concluded that one must have < B > > 200 G (very close to the Hanle 



saturation regime for this line) inside the lanes. This intergranular fields would 
generate most of the Sri depolarization while having no effect in the linear polar- 
ization signals from the molecular lines. It is interesting to point out that when 
the Sri depolarization levels are com puted with the MURaM SSD simulations, 
IShchukina and Truiillo Buenol ( 201lh find that the depolarization levels are far 



too low. This is not surprising as these SSD simulations reach fields of the or- 
der of 20-30 G, very far from the typical 130 G that is needed. The only way in 
which an agreement with the observed depolarization values could be achieved 
was by multiplying everywhere the field strengths in the simulations by a factor 
12. This factor is four times larger than that needed bv iDanilovic et all ( 2010al) 



to match the observed apparent mean Zeeman fluxes of Hinode/SP with those in 
the MURaM simulations. 
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This description of the quiet sun fields as inferred from the Hanle signals is gen- 
erally accepted and no clear alternative exist. However, it is clear th at one would 
like to see it confirm ed by an analysis that is less model dependent fsee lKleint et al.l 
201ltlStenfloll2012l who point out that the analysis made using molecular lines is 
differential in contrast to that made with the atomic lines). The fact that the 
Sri line is so close to the saturation regime has also brought some criticism of 
the a ctual interpretation of the observed depolarization levels (Sanchez Almeida! 
2005). Indeed, very little variation of the Sri polarization signals has been detected 
and this can be used to argue that the tangled hidden field is independent of the 



solar cycle, favoring an or i gin related to SSD action dTruiillo Bueno et aDl2004t 
IVogler and Schussler|[200?l jpietarila Graham et al.ll201ol) . One other explanation, 
of course, could be that the signals in the intergranular lanes change with the 
activity cycle but we do not see the effect because they are always in the Hanle 
saturation regime. 

In any case, the existence of this turbulent unresolved field is well estab- 
lished and one would like to understand if it bears any relation with the inter- 
network fields observed by Hinode/SP and described before. Note that when we 
say unresolved here, we refer to the Hanle observations used in the analysis that 
were obtained over large spatial scales (several arcseconds) and with exposure 
times of the order of one minute or so. Taken at face value, what the Hanle 
measurements need are small scale (arcsecond scale or below) non-vertical field 
patches (the Hanle effect is insensitive to vertical fields), with field strengths be- 
low equipartition with granulation (~ 400 G). It is evident that these are the 
prope rties of the internetwork fields described by Bellot Rubio and Orozco Suared 



( 20121) and lOrozco Suarez and Bellot Rubiol (|201^T and one is tempted to conclude 



that the Hinode/SP fields are al s o the fields corresponding to the Hanle signals. 
iBellot Rubio and Orozco Suare3 (|2012j) specifically mentions this possibility. One 
would be tempted to go further and state that it would be rather strange to have 
the Sun harboring two families of fields with so many things in common but that 
are totally unrelated. Thus, we also favor here the identification of the Hinode/SP 
inte rnetwork fields w ith those that produce the Hanle depolarization signals (see 
also lLites et al.ll2009l) . Perhaps, the only difficulty we encounter with this identi- 
fication is the well known fact that the HIFs (the internetwork transverse fields) 
have a clear preference to be located at the borders of granules (jLites et al.l l2008; 
iDanilovic et al.ll2010bh not in the intergranular lanes as required by the Sri depo- 
larization measurements. However, this could be a minor problem as a granular 
border might be sufficiently close to what is needed. An evident way to test this 
identification would be to invert a volume of the Hinode/SP internetwork observa- 
tions with an inversion code that provides the complete atmo spheric stratification 
such as the SIR code dRuiz Cobo and del Toro Iniestalll992l). Then perf o rm the 
Hanle synthesis as in iTruiillo Bueno et al.l (|2004 and lKleint et al.l (|201ll) for the 



atomic and molecular lines and compare the result with the observations. In this 
exercise some assumptions about the upper layers might be needed to extend the 
retrieved atmospheres over the range of formation of Sri line, but such an extension 
can be reasonably done. 

In any event, it is clear that spatially resolved Hanle depolarization measure- 
ments (|Stenfloll2012l) that tell us where exactly the Sri depolarization signals occur 



in the Sun are urgent. 
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4 Deep magnetograms and 'dead' calm areas: implications 

As commented in Section [TJ dynamo action concentrates in the turbulent down- 
flows as field line stretching and amplification is more efficient there . A dr amatic 
visualization of this can be seen in Figure 1 o f lSchiissler and Vogled (|2008|) . Basi- 



cally all of the intergranular lanes are seen to participate in this dynamo swing. 
One can then expect that the observable effects of an SSD driven by the granu- 
lation would be distributed homogeneously over spatial scales similar to that of 
the granulation after a period of time long compared with the lifetime of an inter- 
granular lane (w hich is similar to that of the granules themselves, 10 minutes, cf. 
iTitle et al.lll989h . In particular, the internetwork fields seen by Hinode/SP must 
be uniformly distributed over granulation scales. Thus, if we are able to somehow 
make a statistics of the location of these fields, the inferred spatial distribution 
must reflect the spatial scale of the granulation and have no voids larger than 
the typical size of one granule (or at least, the probability of occurrence of such 
voids will be small, see below). On the contrary, if one finds areas on the quiet sun 
where many intergranular lanes have existed but none of the expected effects of 
a presumed granular SSD are seen, one can conclude that no such SSD has been 
detected. 

While snapshots from SSDs have been used to make spectral synthesis in- 
cluding the degradation effects of tel escopes and detectors ( Danilovic et al.ll2010al 



IShchukina and Truiillo Buend |201 lh , what one needs is a complete time series 
that allows the study of the generation and disappearance of the dynamo fields 
at current spatial resolutions and polarimetric sensitivities. This study is unfor- 
tunately not available yet. From an observational perspective, we have seen how 
internetwork fields evolve at the solar surface with unprecedented detail. Both 
Hinode/SP and SUNRISE/IMaX have convincingly shown that internetwork fields 
are fed into the surface in the form of emerging small scale loops (rather than in 
the form of, for example, spontaneous appearances of newly created flux patches 
in the intergranular lanes). Hinode/SP with its superior spectral resolution and 



coverage has allowed to study a large number of such events (see ICenteno et al 
20071: IMartmez Gonzalez and Bellot Rubiol 120091; llshikawa et all l2010l: IViticchie 



2012tlGomorv et al-lboToT the latter for ground observations in the infrared) with 



great detail. Basically a horizontal patch is first detected in linear polarization 
that later displays two opposite polarity footpoints that move apart. The dis- 
tance between the footpoints of these small scale loops is typically 1 Mm, the 
lifetime 10 minutes and the magnetic flux of around 10 17 Mx (see th e statistic in 
IMartmez Gonzalez and Bellot Rubiol l2009h . iDanilovic et~aD (|2010bl) . using data 
from SUNRISE/IMaX, identified a large number (thousands) of HIFs that they 
associate with flux emergence in the form of loops. Such a large amount of occur- 
rences emphasize this process as the main source for internetw ork flux at present 
resolu tions and sensitivities. A case in point is that described bv lGuglielmino et al.l 
(2012) who analyzes in detail what is probably the largest quiet sun loop ever ob- 



served. In this case, a maximum footpoint separation of 4.5 Mm is achieved with 
a magnetic flux content of 6 10 17 Mx and a duration of 25 minutes. This quiet 
sun bipole, that h as one order of magnit ude less flux than the smallest ephemeral 
region studied by Ha genaar et al.l (|2003), is arguably not generated by any pro- 
cess that occurs at granular scales. It is known that dynamo simulations generate 
similar loop-like structures (as the horizontal fields described in Section [2] are part 
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Fig. 3 SUNRISE/IMaX deep linear polarization magnetogram integrated over a 30 minute 
period (top). SUNRISE/IMaX deep circular polarization magnetogram integrated over the 
same period (bottom). Both are computed according to the definitions in Eq. 1. 
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of them). But a comparative study of the maximum footpoint separation, flux 
content, etc., is missing. 

A recent study about loop emergence in the quiet sun has recently been pub- 
lished that is relevant for our discussion. Using the two time series of t h e first 
flight day obtained with the IMaX instrument, iMartmez Gonzalez et all (2012) 



studied the emergence of 497 magnetic loops identified in them. They estimate 
an event rate of 0.25 loop h _1 arcsec -2 . If we associate a typical linear size of 2 
arcsec for a granule and a lifetime of 10 minutes, this rate can be translated into 
0.17 loops per granule (it takes 6 granules to get one quiet sun loop). Each one of 
these time series lasted for about 30 minutes, so granulation was efficiently created 
and destroyed over their timespan. Such a large number of detected loops allowed 
them to study the spatial distribution of these magnetic flux emergence processes. 
The result they obtained was that the spatial distribution of loop events was far 
from homogeneous at granular scales. They found what they termed "dead calm" 
areas where simply no loop was seen to emerge during the time series. One could 
argue that these calm areas can be created by chance and that their existence is 
simply a mere coincidence. However, the authors perform a statistical study of the 
likelihood of such voids given their size (70-100 arcsec 2 ) under the assumption of 
a spatially uniform loop emergence probability. They modeled the probability of 
finding one such large circular void with a resulting estimate of 3xl0 -4 . Two such 
dead calm areas were cleanly identified. 

Let us show from another perspective how unlikely this result is. In voids of 
this size, one can fit around 20 typical granules at any given time. As the time 
series covered 3 granular lifetimes, 60 granules existed inside them which would 
have given rise to at least 10 loops at the above rate of creation, but none was 
found. And it occurred in two unrelated regions. How is it possible that if magnetic 
loop emergence is the observable outcome of SSD action, one encounters regions 
where this is not activated? While we find this result highly incompatible with the 
existence of a granular ly driven SSD, we need a solid comparison with simulations 
including all of the possible observational biases to provide a firm answer. 

We must note that a non-uniform distribut ion of quiet sun fields has been 



known for some time. They were identified by ( Lites et al.l I2008L and references 
therein), who pointed out the existence of mesogranular scale voids in the Hin- 
ode/SP map. The preference of internetwork fields to be located at mesogranular 
scales has been clearly demonst rated using the same SUNRI SE/IMaX data as that 
used to identify the voids (see lYelles Chaouche et al.l 1201 lh , but note that these 
voids were much larger tha n a mesogranule and have a scale closer to that of a 
small supergranule (see, e.g. iMeunier et al.ll2007l who give radius in the range of 
8 to 30 arcsec). 

One could argue that whereas no loop emerged in these dead calm areas, 
they were probably not devoid of some subtle form of internetwork field presence. 
However, inspection of the IMaX data with a scaling close to the noise levels 
readily shows that the locations of these voids clearly harbored less activity than 
the rest of the observed area. To prove this point more clearly, we produced two 
time averages of the IMaX time series that are shown in Fig. [3] The two quantities 
that are displayed correspond to deep magnetograms computed as: 

i=l i=l 
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where ey and e^p are quantities inferred from the data that allow convenient 
reduction of the noise when doing the time averages of the otherwise positively 
defined quantit ies. The deep magnetograms in Fig. [3] evidence the same two voids 
as detected by iMartmez Gonzalez et al.l ( 20121) . They are centred at coordinates 
[15,25] and [35,37]. The voids are visible in both the linear and the circular po- 
larization deep magnetograms. Note that the scaling of both magnetograms satu- 
rat es at 3 x 10~ 4 for LP and a t 10 -3 for V . The calibration constants published 
bv lMartrnez Pillet et al.l ( 201ll) would have translated these values into 45 and 5 
G respectively. Measurements below these fluxes are at the limit of state of the 
art imaging magnetographs. These deep magnetograms show that they had re- 
duced levels of magnetic activity and it is correct to refer to them as magnetically 
calm (maybe, not dead). The question, of course, is what was special about the 
many intergranular lanes that populated these regions that prevented them from 
displaying the magnetic activity levels seen elsewhere. 



5 Where do we go from here?: higher sensitivities and higher latitudes 

The si mulations described in Section[T]have reached Re m g [5000,8000] (jPietarila Graham et al.1 
2009b]). However, they fail to provide the magnetization levels needed to explain, 

both, the Hanle ; depolarization signals of the Sri line and the Hinode/SP fluxes (see 

Danilovic et al.l2010al [Shchukin a and Truiillo Buenoll201ll : lOrozco Suarez and Bellot Rubiol 



2012ft . This is often explained arguing that the simulations are still far from the 



magnetic Reynolds numbers of the Sun. But the values that are achieved by them 
are only about one order of magnitude below the expected values at the solar 
surface. In contrast, they reach, at best, similar Re values, which are 10 times 
smaller than what we encounter on the Sun (see Fig. [T] where the values achieved 
by the simulations are marked by an X near the surface) . There are about 5 orders 
of magnitudes in the inertial range below the resistive scale, l-q, populated with 
cells that generate an enhanced turbulent diffusion that are not present in the sim- 
ulations. As mentioned in Section fU recent numerical simulations show that this 
might not be a problem to obtain dynamo action in the low Pm regime. But they 
also tell us that the value of the critical magnetic Reynolds number needed Re„ 
to sustain dynamo action increases sharply when P m < 1. A factor 3-7 increase 
in the value of Re^ is expected. As in the MURaM simulations this number is 
~ 2000, we expect this code to display a dynamo only when Re m > 6000 or larger 
as soon as they use Prandtl numbers in the right ballpark of the problem. Let us 
see the implications for some of the inferences that are made using the available 
numerical simulations. In the case of the MURaM runs, iDanilovic et al. I (|2010al) 



1 1 2 

found a plausible scaling of the saturated field strengths with ~ Re r „ . Their run 
G with Re m ~ 5200 has a mean field of 30 G (at logr w —1). If we scale it 
to a solar value of Re m ~ 10 5 following this square root scaling, we obtain 130 
G. This is the value needed to explain the Hanle depolarization measurements of 
the Sri line (and the most probable f ield strength found in the internetwork by 



lOrozco Suarez and Bellot Rubiol 20121 ) . This nice agreement was already pointed 



out by Shchukina and Truiillo Bueno ( 201ll) . However, for the reasons explained 



above, this is probably a mere coincidence. If this simulation would have been done 
with the same Re m but with P m of, say 0.1 (as in the simulations proving dynamo 
action at low Prandtl numbers), no dynamo action would have been found. The 
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field strength to introduce in the above scaling would have been G instead of 30 
G. 

In discussing the accepted view about the existence of dynamo action at the 
solar surface, we have gone a step further and suggested that the observed spatial 
distribution of quiet sun fields seems to be at odds with a granular ly driven SSD. 
The argument used to make this claim was that in such a dynamo, all the magnetic 
byproducts must necessarily have a uniform spatial distribution at scales above 
that of a granule. We have found, however, that there are voids of magnetic activity 
(dead calm areas) in, both, the average apparent longitudinal and transverse fluxes 
of deep magnetograms and in the distributio n of emerging loops. In fact, we have 
translated the loop emergence rate found bv lMartmez Gonzalez et alj ( 20121) into 
a rate of 1 loop per 6 granules which can also be thought of as 1 emerging loop per 
mesogranule. It is clear that mesogranulation scales are very relevant for the quiet 
sun fields. SSD simulations including them are needed to see if they explain the 
presence of these voids in magnetic activity and emerging loop frequency. Given the 
fact that as we go deeper into the Sun, one reaches higher Re m and P m values, 
it is very tempting to suggest that SSDs acting in a range of convective scales 
larger than the granular ones are those tha t give rise to the presently observed 
internetwork fields. Note that in his work, ICattaneol (|l999f ) already mentioned 
larger scales, such as those of the supergranulation, as possible places where to 
host dynamo activity. This might still be compatible with an SSD at granular 
scales that generates fields much weaker and that have not yet been observed by 
any of our currently available diagnostic techniques. 

We have also suggested that the internetwork fields revealed by Hinode/SP 
and those that generate the Hanle depolarization of some atomi c lines migh t have 
a lot more in common than previous ly thought (see however iLites et al l 2009; 
iBellot Rubio and Oroz co Suare The reason for this identification is based 

on the results from lOrozco Suarez and Bellot Rubiol ( 2012h and lBellot Rubio and Orozco Suared 
(|2012h who have shown that the quiet sun fields share the same field strengths, 
inclinations and azimuths than those needed by the fields detected with the Hanle 
effect. This is, of course, compatible with a continuous spectrum of fields in which 
the ranges corresponding to the two types of fields (the internetwork and the 'hid- 
den' Hanle fields) simply overlap over a much larger fraction than thought so far. 
Of course, the internetwork fields will also include a fraction of vertical kG fields 
that contributes nothing to the Hanle depolarization and the weak granular fields 
that depolarize the light in molecular lines never make an imprint in the Fel Zee- 
man lines observed by Hinode/SP at current sensitivities. But, of this continuous 
distribution of fields, the range of hG strengths with largely transverse orienta- 
tions and spatially organized at granular scales, contributes simultaneously to, 
both, the Zemman and the Hanle observations. If this result is confirmed (as, e.g., 
with the SIR inversions and the Hanle computations mentioned in Section [3]) an 
important step to clarify the currently complex debate of the nature of the quiet 
sun magnetism would be achieved. We also want to stress that the evidence that 
the internetwork field component is composed of largely transverse fields renders 
the debate about the mean value of < \B C ^ P \ > (or of < \Bl\ > for that matter) 
obsolete. Stokes V (longitudinal) magnetograms of the internetwork simply show 
a rather incomplete picture of these fields. 

Note that many of the results commented above have benefited from high po- 
larimetric sensitivities. Hanle measurements have always been very demanding in 
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polarimetric accuracy. The deep ma gnetograms of SUNRISE/IMaX and th e long 
integrations with Hinode/SP used bv lOrozco Suarez and Bellot Rubiol (j2012ri both 
were at the 1CP 4 polarimetric sensitivity. This is not a coincidence. Much of the 
future progress will be achieved with sensitivities in this range. Those more regu- 
larly reached in present day observations, 10 -3 , are due to instrumental limitations 
that have nothing to do with physical processes in the Sun. Polarimeters observ- 
ing the solar photosphere with sensitivities of 10~ 4 and sub-arcsecond resolutions 
using the Zeeman and the Hanle effects will consolidate (or refute) many of the 
aspects commented here. The need for high spati al resolution observations of the 
Sri depolarization cannot be emphasized enough (jStenfloll2012f ). These targets de- 



mand large apertures similar to those plan ned for future facilities such as the 4m 
class ground-based telescopes (ATST, EST lKeil et~al1l201lt ICollados et al.ll2010l) 
and the Japanese led Solar-C mission (1.5m aperture). All these facilities will 
likely have to be used outside of the diffraction limit t o pursue high sensitivity 
spectropolarimetry thanks to their large collecting areas iKellerl (|l999l) . 



Finally, observing regions of the Sun hardly reachable from the ecliptic will also 
help to clarify the nature o f the quiet sun magnetism. The ESA-led Solar Orbiter 
mission (jMiiller et al.ll2012l) will carry on-board a ma gnetograph (the Polarimetric 
and Helioseismic Ima ger, PHI; iGandorfer et al.|[2oTlh similar to that of SDO/HMI 
( Scherrer et alJ 120121) and will observe the dynamics of the solar poles from an 



inclination of 35° with respect to the solar equator. Observing the poles is crucial 
in this discussion because they represent the regions at farther distances from 
the activity belts in the Sun. In the absence of an SSD mechanism working at 
the solar surface, the origin of the internetwork fields can only be explained as 
a result of the cascading d own towards the smallest scales (see the discussion in 
ISchiissler and Vogrerl l2008) of the global dyn amo fields. This effect is inevitably 
present on the Sun, but whether it affects only to network fields or to internetwork 
ones can best be discerned by observing their latitudinal properties with good 
spatial resolution and sensitivity. 
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